Adatoms in graphene 



A. H. Castro Neto, Valeri N. Kotov, and Vitor M. Pereira 
Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA 

Johan Nilsson 

Instituut-Lorentz, Universiteit Leiden, P. O. Box 9506, 2300 RA Leiden, The Netherlands 



oo 

o 

O 

(N 

O 

Q 



CZ3 



N. M. R. Peres 

Centra de Fisica e Departamento de Fisica, Universidade do Minho, P-4710-057, Braga, Portugal 

Bruno Uchoa 

Department of Physics, University of Illinois at Urbana- Champaign, 1110 W. Green St., Urbana, IL 61801-3080 

(Dated: December 11, 2008) 

We review the problem of adatoms in graphene under two complementary points of view, scat- 
tering theory and strong correlations. We show that in both cases impurity atoms on the graphene 
surface present effects that are absent in the physics of impurities in ordinary metals. We discuss 
how to observe these unusual effects with standard experimental probes such as scanning tunneling 
microscopes, and spin susceptibility. 
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It has been known since the discovery of graphene [1, 2] 
that disorder plays a fundamental role in determining its 
physical properties [3]. It was understood early on that 
the linear dependence of the graphene conductivity with 
electronic density, or gate bias, at large densities requires 
the presence of long range scattering, either in the form of 
Coulomb interactions [4-6] , or correlated ripples [7] . Ex- 
periments with K atoms in ultra high vacuum conditions 
have shown that the conductivity is affected by ionized 
impurities [8], although that might depend strongly on 
the screening properties of the environment [9]. In fact, 
because of its electrical sensitivity to adatoms, it is pos- 
sible to use graphene as a single molecule detector [10]. 

Although unwanted adatoms can be a nuisance for 
transport properties, one can also think of them as a 
way to tailor graphene, in order to create new many- 
body states that do not show up in pure graphene [11]. 
The situation here is akin to what is found in ordinary 
metals which, when in pure form, are weakly interact- 
ing and inert, but show extraordinary properties such as 
magnetism, Kondo effect and superconductivity, when 
doped with transition metal atoms with strongly inter- 
acting inner electronic shells [12]. 

The advantage of graphene in comparison to ordinary 
three dimensional (3D) metals (where adatoms are intro- 
duced by alloying, which is generically a random process) 
is that graphene is an open surface, and hence adatoms 
can be manipulated with the use of atomic force mi- 
croscopy to obtain structures with atomic precision [13]. 
Only recently, has scanning tunneling microscopy (STM) 
been used to study adatoms in graphene [14], and still 
very little is known about the spectroscopic properties of 
these adatoms. 

Historically there have been two main approaches for 



the understanding of impurity atoms in ordinary met- 
als. The most traditional one, championed by J. Friedel, 
follows the arguments of traditional scattering theory in 
terms of the scattering phase shifts of the electrons in the 
presence of an external potential [15]. In this kind of ap- 
proach, the electrons of the impurity atom and the metal 
host are treated at equal footing and the many-body as- 
pects of the problem enter only through Pauli's exclusion. 
As shown by Anderson [16], this kind of approach is in- 
adequate for atoms with inner shell electrons, where the 
Coulomb interaction is large and has to be taken into 
account from the outset. In this case, one distinguishes 
between the host and the adatom electrons. These two 
approaches are essentially complementary and serve to 
address different issues in the problem of impurity atoms 
in metallic hosts. 

We will discuss these two approaches in the context of 
adatoms in graphene. In the fist part of this paper we 
discuss the Friedel approach. We are going to show that 
the Coulomb problem in graphene has features similar to 
the ones found in the 3D relativistic Coulomb problem. 
In particular, the famous "fall to the center effect" [17] 
that would occur in atoms with atomic number Z greater 
than l/oc, where a — e 2 /{he) is the fine structure con- 
stant, can occur in graphene even in the presence of a 
single proton if the dielectric constant eo is of the order 
of 1. This effect can be observed by STM spectroscopy 
by studying the local density of states (LDOS) close to 
an impurity. It is interesting that graphene allows for 
an opportunity to study an effect that was predicted to 
occur only in heavy ion collisions, or in the event horizon 
of a black hole [18]. 

In the second section, we present the generalization of 
the Anderson impurity model to graphene, taking into 
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account the local Coulomb repulsion (Hubbard term) of 
inner shells. We show that the Anderson problem in 
graphene has features that are rather different from the 
ones found in ordinary metals [16]. In particular, we ar- 
gue that transition metals that would not be magnetic in 
ordinary metals can become magnetic in graphene due to 
its peculiar properties. In the last section we include our 
conclusions. 



FRIEDEL'S WAY [19-22] 

In what follows, we use the effective description of the 
electrons at low energies and long wavelengths [23] by ex- 
panding the electronic dispersion around the Dirac point 
in the Brillouin zone (K and K' points) and obtaining the 
two-dimensional Dirac (we use units such that H = 1): 



H = is F & p 



(1) 



where v-p w 6 eV A is the Fermi-Dirac velocity [23], 
o i [i = 1,2,3) are Pauli matrices, and pi (i = x,y) is 
the two-dimensional momentum operator. Consider the 
problem of a single Coulomb impurity with charge Ze 
sitting at the origin. The wavefunction of the problem 
obeys the equation: 



up {-ia ■ V - -)#(r) = E *(r), 



(2) 



where &(r) 
and 



(ipA( r ),4'B( r )) is a spinor wavefunction, 
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Ze 2 



(3) 



is the dimensionless coupling constant of the problem. 
Hence, the coupling constant depends directly on the di- 
electric properties of the medium where graphene is sit- 
ting [24]. 

Notice that Eq.(2) only involves one of the Dirac 
cones. This assumption is justified on the basis that the 
Coulomb interaction is long-ranged, and hence singular 
in momentum space (it behaves as 1/p as p — > 0). There- 
fore, scattering is strongly localized around one of the 
cones. We have checked this assumption numerically by 
solving the Coulomb problem in the tight-binding model 
in the lattice [19]. 

Observe that the total angular momentum, J z = L z + 
cr 2 /2 [25], is conserved and we can therefore diagonalizc 
the problem in the basis of momentum states: 



1 /jU-itecpffr) 



(4) 



in cylindrical coordinates, Eq. 
±1/2, ±3/2,...) 

( £ + g/r -(dr + j/r) 
\(9 r -j/r) e + g/r 



(2) becomes: (j 



(5) 



The solutions for this problem can be written as: 

^'( r ) = X] C n U nfn(r) , li±l = 
n=±l 

where s x = sign(a;), 



l_f ^\ j±X\ 
m\s gj ^/\jTX\J , 
(6) 



and f\(r) is the solution of the equation: 
2 , 2ge A(ATl) 



(7) 



d z r f±i{r) 



/ ±1 (r)=0, (8) 



which is the radial equation for the 3D Coulomb problem 
[17] with the energy, E, replaced by e 2 , and the angular 
momentum quantum number, t, replaced by A. 

We immediately notice from (7) that A can be a purely 
imaginary number when g > j. The presence of "imag- 
inary angular momentum" channels in the problem can 
be initially troublesome but they are mathematically al- 
lowed. Notice that, since j is a half-integer, the smallest 
value of g for this to happen occurs for g = g c = 1/2. 




FIG. 1: (color online) Phase shift of the Coulomb problem in 
the undercritical regime as a function of coupling g. 

For g < g c , the so-called undercritical regime, eq. (8) 
can be solved as (see (6)): 

<Pj(r)/J\fj =u + Fx-i(-g,p) + s ge u-Fx(-g,p), (9) 

where A/} is a normalization factor, F^(j], p) is a Coulomb 
wave functions [26, 27], p = \e\r, and g — s e g. At 
long distances, the asymptotic behavior of this function 
is given by [17]: 



\(-g,p) ~ sm(p + glog(2p)+ti x (g) 



(10) 



3 



where i? A (<?) = -Af + arg[r(l + A - ig)] . Hence, in (9) 
we have asymptotically: 



sm 



p + <?log(2p) + arg(u+e 



«X-1 



(11) 

One notices the logarithmic phase shift which is tradi- 
tional in Coulomb problems. Observe that the phase 
shift is energy independent in this case, but it is strongly 
dependent on the value of the coupling g. This result has 
implications to the so-called Friedel sum rule. In Fig.l 
we plot the value of the phase shift as a function g. 

The presence of the adatom modifies the local density 
of states (LDOS): 

N(e,r) = 5> £ (r)| 2 <5( e -i?) 

E 

oo 



j=-oo 



where 



nj {e,r) = r-^fWI 2 



1^ 



-l,,^ (r) |2 



2gF A F A _ 1 /|j|] (13) 



is the contribution of each momentum channel to the 
LDOS. In the limit g ^ the Coulomb wave func- 
tions reduce to Bessel functions [26] , and one obtains the 
free DOS [23] if the normalization is chosen as Mj 2 = 
2tt 2 X 2 /j 2 . It is clear from (13) that particle-hole symme- 
try is broken, and as e — > we have 



N(e,r) 



2A 



(14) 



and A < 1/2 for \j\ = 1/2. At higher energies, or far away 
from the adatom, we obtain the bare density of states of 
graphene, namely, N (e,r) ~ |e|. 

For g > g c , the so-called supercritical regime, the solu- 
tion of (9) is of the form, ifij(r) = C\ (pip{f)+C2 ^P-if3{r) , 
where C\ t 2 are to be set by the boundary conditions at 
short distances and 

<Pip(r) = u+F l/ 3_ 1 (-g, p) + s jgt :U^F tl3 (-g, p), (15) 



(3 = —i\ = \J g 2 — j 2 for those j's such that \g\ > 
and 



u±i 




SgVj T W 



(16) 



The issue of boundary conditions at short distances is 
quite interesting in this case. Notice that as r — * 0, the 
behavior of the Coulomb wavefunctions with imaginary 
index is rather peculiar: 



F ip ^(-gt,p^0) ~ ~ oos(/31n(p)) , 



(17) 




FIG. 2: (color online) Phase shift of the Coulomb problem as 
a function of energy E in the supercritical regime (g > g c ). 



and, hence, the wavefunction oscillates endlessly, with 
an infinite number of zeros, as one approaches the origin. 
This is a rather puzzlingly feature since in basic quantum 
mechanics the energy of a state is directly related to the 
number of zeros of the wavefunction. 

This effect is known as the "fall of a particle to the 
center" [17] and occurs when the external potential is too 
strong [18]. Notice that the condition g > g c = 1/2 is 
equivalent, in terms of (3), to: 



Z > 



1 



where 



a G 



2a ( 



(18) 



(19) 



is graphene's "fine structure constant" [28, 29]. In rel- 
ativistic quantum mechanics this effect shows up as the 
"fall" of the electron into the nucleus when the atomic 
number is too high. As the electron falls, a positron 
is emitted leading to the charging of the vacuum. The 
analogy in graphene is clear: if the potential is attractive 
to the electrons, the electrons become strongly bound to 
the adatom and a "hole" is ejected. Unlike the problem 
in free space, the problem at hand occurs in a lattice 
which always has a natural ultraviolet cut-off, namely, 
the lattice spacing or in this case the distance of the 
adatom from the carbon atoms, ciq. Hence, this indef- 
inite oscillation of the wavefunction is controlled by a 
short distance cut-off do . The vanishing of the wavefunc- 
tion at r < clq translates [30] into the boundary condi- 
tion: ipf(ao) = ipf(a ). This allows us to fix the value 
of Ci / C2 and determine the total scattering phase shift 
in the supercritical regime. Unlike the undercritical case, 
the phase shift is a strong function of the energy, as shown 
in Fig. 2. 
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(a) 




(c) 



FIG. 3: (color online)LDOS as a function of energy at differ- 
ent distances r from the adatom: (a) r — ao; (b) r = 2oo; (c) 
r = 3oo; (d) r — » oo. 



We can now proceed and compute the contribution of 
the overcritical j's to the LDOS: 



n 3 {e,r) 



1 g>) + s gj Re[e i! %f(p)] 
2n 2 r q!(oc) + s ej Re [e i25 ^j i '(oo)] ' 



(20a) 



\Fi P \ 2 + \F^\ 2 + 2\j\ ReiF^F^^/g, (20b) 

(20c) 



where, 

eh 

g 1 / ee 2F tf3 F lf ^ 1 + + F%_ x )/g. 

Eqs. (13) and (20) determine the LDOS for any coupling 
strength, g, which can be summarized as 

N(e,r)= Y, »i(e,r)+ £ n ^ r )- ( 21 ) 
b'|<|s| |j|>lffl 

In Fig. 3 we show the LDOS as a function of energy at 
different distances from the adatom. One can clearly see 
formation of resonances in the LDOS that are directly 
associated with the positron creation of the actual rela- 
tivistic problem. The possibility of testing ideas of rel- 
ativistic physics in a carbon-based material is unique to 
graphene due to the robustness of its Dirac spectrum. 

In the previous discussion we have completely disre- 
garded the Coulomb interaction between the electrons. 
There are three types of electron-electron interactions 
one should worry about: electron-electron interactions 
within graphene electrons, interactions between graphene 
electrons and adatom electrons, and between electrons 
in the adatom. Electron-electron interactions within 
graphene, and between electrons in graphene and in the 
adatom, can change some of the quantitative features dis- 
cussed here, as we shall see below. The electron-electron 
interactions within the adatom can lead to new physics, 
namely, the possibility of creating a local magnetic mo- 
ment in the adatom. This kind of situation is not con- 
templated in Friedel's way and will be the subject of next 
section. 



Let us consider the problem of screening in graphene. 
Within the random phase approximation (RPA), in the 
static case, the screened potential of an external impurity 
is given by [31]: 



V s {q) = 



V (q) 



l-U(0)( q )U(q) 



(22) 



where Vo(q) = ZU{q) is the Fourier transform of the 
impurity Coulomb potential in 2D, U(q) — 2ire 2 jq the 
Fourier transform of the bare electron-electron interac- 
tion, and IT )(g) is the static polarization function. In 
the above expression we have disregarded the modifica- 
tion of the dielectric function due to electron-electron 
interactions [32]. In a metal, we have n(°)(<7 -> 0) oc 
—N(Ep), the density of states in the Fermi level, and 
hence V s (q — » 0) — » constant, indicating that the poten- 
tial is completely screened. 

Another way to understand screening is by comput- 
ing the total induced charge density. The local particle 
density in RPA is given by [31]: 



Sn(r) 



Il( \q)Vo(q) 
aq- z^ttx ; r~rr , e 



l i-n(»)(?)%) 

and hence the total induced particle density is: 



(23) 



SN = J drSn(r) 

n(°>( s - 



0)V o (q^0) 



l-n(°)(t)f->0)Z7(g->0) 



(24) 



and one immediately sees, in very general terms, that if 
Vo(q — > 0) and U(q — > 0) are divergent and II' ' (q — > 0) is 
regular, as in the case of a normal metal, the RPA result 
is that 



SN, 



metal 



(25) 



that is, the local impurity is completely screened. 

Let us apply the same reasoning to graphene. Because 
RPA is valid in linear response, it only applies to weak 
coupling, that is, in the undercritical regime, g < g c . 
The polarization function in graphene has been calcu- 
lated many times before [33-35] , but the argument above 
only requires knowledge of the polarization function at 
long wavelengths. The asymptotic behavior of the polar- 
ization function is: 



n (o) (g^0,/i = 0) 



q 

4z/ F 



(26) 



for ii = (that is, when the Fermi energy is at the Dirac 
point), and 



n(%^o,^o) 



(27) 
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when /i^fl. Replacing these results in (22) and (24) 
give: 



2ne 2 



1 + §a G q 



(28) 



for jU = 0, and 



2^Ze 2 
~ — r^- 

9 + 9o(m) 



(29) 



for fi ^ 0, where <jo = 4ac|/i| is the Thomas- Fermi 
screening length. This result indicates that graphene 
screens like a metal whenever \x is finite but when the 
Fermi energy is at the Dirac point there is a residual 
unscreened charge, 



SN uns . = 
Z* = 



Z - 5N(n = 0) 
Z 



1 + 5«G 



(30) 



in accordance with (28) showing the Coulomb potential 
remains the same even with screening and that the only 
modification is the value of the local charge that changes 
from Z to Z*. Notice that the limits of \i — ► and q — > 
do not commute in graphene. In fact, the limit of /x = 
and q — > is rather pathological. 

In the overcritical regime the problem of electron- 
electron interactions becomes of fundamental impor- 
tance. As we have seen previously, the "fall of the parti- 
cle to the center" indicates that large amounts of charge 
accumulate in the vicinity of the impurity. This accu- 
mulation of charge implies strong electron-electron in- 
teractions beyond the description of RPA. Instead, the 
screening becomes non-linear [36, 37] and the Coulomb 
law is modified. This regime, however, is highly non- 
perturbative in nature and it is not clear that the ap- 
proximations that are usually used in this regime, such 
as the non-linear Thomas-Fermi theory, are really appli- 
cable. 



ANDERSON'S WAY [38] 

In Anderson's approach we distinguish between the 
graphene electrons, described by Hamiltonian (1), and 
the electrons localized in the inner shell of the adatom 
which have an Hamiltonian: 



(31) 



where f a (/t) annihilates (creates) and electron with spin 
a =T,I, £q is the bare energy of the local level, and U 



is the energy cost to doubly occupy that level. The in- 
teraction between graphene electrons and local level is 
described by a hybridization term: 



H v = Vj2[f^BA0) + H.c] 



(32) 



where V is the hybridization energy, and l 4>B,a{ r ) is one 
of the components of the spinor defined earlier. In (32) 
we have assumed that the adatom hybridizes on top of a 
carbon atom (in this case, one that is sitting on sublatticc 
B). This assumption can be easily relaxed and the only 
difference is a form factor that takes into account the 
geometry of the position of the adatom. Since we will 
be interested in s-wave scattering only, the details about 
this form factor will be irrelevant in what follows. 

The complexity of the problem comes from the U term 
that describes double occupancy of the adatom. In the 
absence of this term the problem is quadratic and can be 
diagonalized exactly [39]. From the point of view of the 
adatom, graphene behaves as a heat bath that damps the 
electronic motion in the adatom. As a result, the local- 
ized electrons acquire a self-energy, T,ff(ui), as a func- 
tion of frequency ui. The real part of the self-energy 
is associated with the shift of the local level energy, 
e/ = e — 5i£//(eo), while the imaginary part, 3S//(o>), 
gives the decay rate, or broadening, of the adatom level. 
A straightforward calculation gives: 



A 2 



A 2 



S&ffiw) = -V 2 ^6(k-\u\). 



(33) 



where A (w vp /a ) is a ultraviolet cut-off. It is imme- 
diately clear that the problem at hand is rather different 
from the case of adatoms in ordinary metals. For one, in 
ordinary metals the damping of the local level is essen- 
tially energy independent and proportional to the DOS 
at the Fermi level of the metal. In graphene, the DOS 
vanishes linearly with energy and this is immediately re- 
flected in (33) where the imaginary part of the self-energy 
behaves as \u>\ and hence vanishes at the Dirac point (by 
Kramers-Kronig the real part necessarily has to behave 
as uo ln(w) as uo — > 0). This unusual broadening has direct 
consequences for the LDOS. 

Notice that the LDOS of the adatom (at the posi- 
tion of the adatom) changes from a Dirac delta function, 
pfJ t<T (u) = S(uj - e ) to: 



PffA u ) 



where 



A\w\Z(u) 



7T [w - e a Z{uj)} 2 + (A|w|Z(w) + 0+) 2 



A = 7T 



A 2 



(34) 



(35) 
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is the dimensionless hybridization parameter, and 

Z~ 1 (uj) = 1 + (V 2 /A 2 ) In (|A 2 - lo 2 \/uj 2 ) (36) 

is the quasiparticle residue of the local electrons. One of 
the most striking features of (34) is that, although it is 
localized at the energy of the level, it has a strong non- 
lorcntzian shape with a substantial tail, something that 
should be easily measured by STM. Moreover, the quasi- 
particle residue vanishes at the Dirac point indicating 
that the local electron is not a well defined excitation at 
the Dirac point, being completely merged with the other 
graphene electrons. These are features that do not show 
up in ordinary metals. 

The simplest way to include the Coulomb U into 
the problem is, following Anderson's original paper [16], 
to perform a mean-field factorization of the interaction 
term, namely, n^ni — > (n-f)nj +n^{n{), where 



«) = r dupff^u) (37) 

J-A 



(n, 



is the average adatom occupation that has to be calcu- 
lated self-consistently. At the mean field level, the only 
effect of the interaction is the shift the level energy in a 
spin-dependent way: 



U(r 



(38) 



and hence, in order to calculate (37) we can use the non- 
interacting result (34) by simply replacing the bare level 
energy, eo, by is renormalized level (38). The integration 
in (37) is straightforward and one obtains (Z(e a ) = Z a ): 



n„ 



(Z„ 2 + A 2 ) 



, /i . 1 

(— e a ) + -arctan 




where 
Qcr = Z a — In 

7T 



w a {E a y< 



where 7 = sign(^), and 

E a = \J (e a - /iZ, 
W a = J{AZ^ + 



1 / AA 

arctan —, 

tt \AZ a l +e a 



-^2 



u 2 A 2 



A 2 A 2 . 



(40) 

(41) 
(42) 



The solution of the coupled set of equations (37), (38), 
and (39) gives the spin population of the local level in 
the adatom at the mean-field level. Notice that the fact 
that the energy depends on the spin orientation through 
(38) allows for states with different spin occupation and 
hence with a local moment, (m) given by: 



(to) = (n T ) - (n x ) 



(43) 



and the states of the adatom can classified as either mag- 
netic, (to) 7^ 0, or non-magnetic, (to) = 0. One would 
like to know the boundary line between these two states. 



Following Anderson it is convenient to define two di- 
mensionless parameters, x = AA/U and y = (p — eo)/U, 
that measure the degree of hybridization between adatom 
and graphene, and the distance in energy of the level 
from the Fermi energy, respectively. In Fig. 4 we show 
the boundary line between magnetic and non-magnetic 
states as a function of the parameters of the problem. 

Fig. 4 shows that, just like in a normal metal, one 
needs sufficiently strong Coulomb interactions to gener- 
ate a magnetic state. Nevertheless there are substantial 
differences between this diagram and the one obtained 
in the problem of impurities in ordinary metals [16]. In 
an ordinary metal the whole boundary line is essentially 
scale invariant and depends only on the values of x and y. 
In the case of graphene this is not so and the boundary 
line also depend on the cut-off A indicating the sensitiv- 
ity of the problem to the details in the ultraviolet. This 
result could be advanced from the form of the self-energy 
(33) since it is explicitly dependent on A. Therefore, it 
would be important to investigate the problem again with 
the full band of graphene, not only the linearized effec- 
tive theory (1). Moreover, clearly the adatom breaks the 
particle- hole symmetry around the y = 1/2 line. This 
comes about because the energy of the localized level 
and the Fermi energy can be either above or below the 
Dirac point. On the one hand, then eo > we see that 
there is a magnetic state even when y < 0, that is, when 
6q > /1 and the level should be nominally empty. On the 
other hand, when eo < there is magnetic state when 
y > 1, that is, when \i > e + U and the level should 
be nominally doubly occupied. These results contradict 
our intuition based on the physics of normal metals but 
once again this can be traced back to the form of the 
self-energy (33) which is non-lorentzian with long tails 
at high energies. These tails allow for formation of local 
moment even when the position of the energy level is far 
away from the chemical potential. As a result, it is eas- 
ier to generate magnetic moments in graphene than in 
an ordinary metal and one would expect that atoms that 
are not naturally magnetic on a metal host to become 
magnetic in graphene. 

The application of a potential V g through an electric 
field via a back gate [2] shifts the chemical potential /1 and 
moves the magnetic state of the impurity in the vertical 
direction, that is, the value of y in Fig. 4 for a fixed 
value of x = AA/U. In this way, the magnetization of the 
impurity can in principle be turned on and off, depending 
only on the gate voltage applied to graphene. 

The impurity susceptibility is defined by \ — Mb(«t — 
ni)/B, where hb is the Bohr magneton, and B is the 
applied magnetic field. In the presence of a field, the 
energy of the impurity spin states changes from (38) to 



e a = Cf — ajjLsB + Un- a . 



(44) 



In the zero field limit, the magnetic susceptibility of the 
impurity \ = X) CT a (d«o-/di3) B=0 can be calculated 
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FIG. 4: Boundary between magnetic and non-magnetic impu- 
rity states in the scaling variables x and y for eo > (a) and 
e < (b). Circles: [eo|/A = 0.029, V/A = 0.14; Squares: 
e /A = 0.043 and V/A = 0.14; Triangles: | e 1 /A = 0.029, 
V/A = 0.03. The upturn close to y = 1 and x — > on panel 
b) is not visible in this scale when V is very small (triangles). 
See details in the text. 
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FIG. 5: (color on line) Spin polarization, ri|,and n± and mag- 
netic susceptibility vs. fi for — eo|/A = 0.029 and V/A = 0.14. 
Left panels: x = 11 (dashed curves), and x — 5 (solid). The 
vertical line marks the position of the level, eo = 0.2 eV. On 
the right: Comparison between eo = 0.2 eV (black curves) 
and eo = —0.2 (brown) at x — 0.43. 



straightforwardly from Eq. (39), 



X = 



2 \ dn 

" Mb 2^ ~dT 



1 - u 



dri- 



de- 



i-u 2 



dn- a dn„ 



(45) 



which is shown in Fig. 5, together with the spin occupa- 
tion, as a function of chemical potential. The possibility 
of tuning the magnetism of an adatom with an external 
electric field is unique to graphene. 

The discussion above does not include the Kondo ef- 



fect, namely, the magnetic screening of the adatom mo- 
ment by the graphene electrons [40, 41], and hence it can 
be thought to be valid for temperatures above the Kondo 
temperature, Tk- The situation here is somewhat remi- 
niscent of the case of magnetic impurities in d-wave su- 
perconductors [42] although in graphene the chemical po- 
tential is free to vary with the applied gate voltage while 
in a superconductor the chemical potential is pinned at 
the superconducting gap. Moreover, the quasiparticle ex- 
citations in superconductors have electron-hole character 
and essentially do not carry charge, while in graphene 
they do carry charge. More interesting, perhaps, is the 
fact that the anomalous broadening of the local level 
puts graphene closer to a mixed valence system, that 
is, it allows for different oxidation states. In this case, 
charge fluctuations, which are not usually treated cor- 
rectly in Kondo- type, spin-only, models, cannot be ne- 
glected and one has to deal with the Anderson Hamilto- 
nian directly. It is clear that the calculations described 
here are only the first step towards the understanding 
of the magnetic properties of atoms in graphene. With 
the use of use of the current STM technology one can 
envision the tantalizing possibility of creating and con- 
trolling, at the atomic level, Kondo lattices in the surface 
of graphene. Since the correlated state of Kondo lat- 
tices, the so-called heavy fermion state, is still one of the 
most intriguing problems in modern condensed matter 
research [12], graphene presents a possibility for answer- 
ing many of the questions that cannot be answered in 
ordinary metal hosts. 



CONCLUSIONS 

We have shown that the scattering analysis of the 
Coulomb problem in graphene has two regimes: under- 
critical when g < g c — 1/2 and supercritical when g > g c . 
In these two regimes the behavior of the scattering phase 
shifts and the wavefunctions are very different. In the 
undercritical regime the physics is more or less conven- 
tional, the Coulomb potential breaks particle-hole sym- 
metry, and leads to an asymmetry in the LDOS close 
to the adatom. In the supercritical regime resonances 
show up in the LDOS, mimicking the "fall to the cen- 
ter" of an electron with concomitant positron production 
and, hence, the charging of the vacuum. It is interesting 
to note that the transition between these two regimes is 
controlled by the charge of the impurity Z and the dielec- 
tric properties of the environment, eo. The supercritical 
regime may be reached by eliminating dielectric material 
such as H2O, and working with alkaline atoms such as 
Ca, that have high oxidation states. 

When electron-electron interactions in the adatom are 
included via the Anderson impurity model, we have 
shown that, at the mean- field level, a picture very dif- 
ferent from ordinary metal hosts emerges. Due to the 
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anomalous broadening of the local level in graphenc, 
magnetic states that are nominally not allowed can be- 
come magnetic. This opens the possibility for the mag- 
netization of adatoms that are not usually magnetic in 
normal metals. Since graphenc is an excellent conductor 
and magnetism can theoretically be achieved, one can 
envision using graphene for spintronics applications [43] 
where graphene works as the conductor of charge, and 
the adatoms work as magnetic memory. In this case, we 
would have the logical and memory materials embedded 
in the same matrix, but more important, and unlike di- 
luted magnetic semiconductors [44], the position of the 
adatoms can be controlled with atomic precision with a 
STM. 

It is clear from the previous discussions that, no mat- 
ter how we look at the Coulomb problem in graphene, 
either from the perspective of scattering theory, or mag- 
netism, it has many unusual features that do not appear 
in conventional metals. Properties such as supercritical 
behavior and "fall to the center" physics, and magne- 
tization in nominally forbidden regions of the magnetic 
diagram, are not intuitive and emerge in ordinary matter 
only when those are subject to extreme conditions such as 
enormous electric fields. Nevertheless, in graphene such 
extraordinary properties can be measured with standard 
condensed matter experimental probes such as scanning 
tunneling microscopy. Moreover, such unusual properties 
may be helpful in order to generate strongly-correlated 
states by chemically modifying graphenc, a field that is 
still in its infancy. 
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